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CN ■ Abstract 

The gluon distribution is dominated by the hard pomeron at small x and all Q , with no soft- 
pomeron contribution. This describes well not only the DGLAP evolution of the hard-pomeron part 
of F%{x,Q ), but also charm photoproduction and electroproduction, and is consistent with what is 
[ known about the longitudinal structure function. 

<N 

There is still no agreed fundamental explanation for HERA's striking discovery, the surprisingly large 
l/") , rise with increasing 1/x of the proton structure function F%{x, Q 2 )- This rise is seen to become more 
marked as Q increases. The conventional viewM 

is that the steeply-rising component is absent 
Tj" | at small Q 2 and as Q 2 increases it is generated through pQCD evolution. We have argued^] that 
this view is mathematically suspect, since it relies on an expansion of the splitting matrix which is 
likely to be unsafe W because it induces singularities that are almost certainly not present in the exact 
matrix. Instead, therefore, we maintain that a rapidly-rising term is present already at small Q 2 , so 
that pQCD evolution does not generate this term, but merely makes it become more prominent as Q 2 
increases. 
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Our view is supported by ZEUS measurements^] of the charm structure function F£(x, Q 2 ), which 
already displays a strong rise with increasing 1/x even at small Q 2 . Figure 1 shows the data at 
Q 2 = 1.8 GeV 2 ; the lines are the conventional fitM to the data (MRST) and our own calculation. 
Indeed, HERA data for photoproduction and electroproduction of charm have the striking property^] 
that at each fixed Q 2 they vary as the same power x~ e °, with eo ~ 0.4. This behaviour is not 
widely appreciated as the data are normally shown on a log-linear plot rather than a log- log plot. By 
definition, we say that it is associated with the exchange of an object known as the hard pomeront'''] . 
A term with the same power is present also in the light-quark contribution to the proton structure 
function. To a good approximation, the hard-pomeron coupling to the charmed quark is found to have 
the same strength as to each of the light quarks. Once this assumption of flavour blindness is made^J 
the hard-pomeron component of the complete structure function F%(x,Q 2 ) immediately provides a 
successful zero-parameter description of its charm component F^x^Q 2 ) at small x: 

F%(x, Q 2 ) = A c (Q 2 ) 1+eo (1 + Q 2 /Qo)" 1_eo/2 *" eo (1) 
with Qq ~ 3 GeV and i c ~6x 10 -4 . We define the charm cross section 

47r 2 a F 



v c (W) = ^^FC(x,Q 2 ] 



x=Q 2 /(W 2 +Q 2 ) 



With (1), 

a c (W) = 0.066 W 2t °/(1 + Q 2 /9.1) 1+e °/ 2 (2) 
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Figure 1: charm structure function dataK at Q 2 = 1.8 GeV 2 with curves from MRSTW andt 8 ! the 
form (1) 



where the units are fib and GeV. This expression corresponds to the thin lines in the plots of figure 2. 
These lines are not a fit to these data; they are calculated from the fit to the hard-pomeron component 
of the complete structure function F2(x,Q ), making the flavour-blindness assumtion. 

The hard-pomeron contribution to the complete F2(x,Q 2 ) is the same, with A c replaced with 
A s=s 1.5 x 10~ 3 . We have shownt^l that DGLAP evolution, with a gluon structure function that 
is dominated at small x by hard-pomeron exchange alone, produces a Q 2 dependence for the hard- 
pomeron part of F2 that agrees numerically with this. Our procedure, which can be applied only 
to the hard-pomeron component, gave almost identical outputs for LO and NLO evolution. A good 
numerical fit to the output of the DGLAP evolution for the small- x behaviour of the gluon structure 
function is 

xg(x,Q 2 ) ~ 0.95(Q 2 ) 1+£o (l + Q 2 /0.5)- 1 - en/2 x- £o (3) 
This fit is valid for Q 2 between 5 and 500 GeV 2 . 

We use the gluon distribution (3) to calculate charm production in leading-order pQCD and compare 
the result with (1). First, we calculate simple photon-gluon fusion in lowest order, for which the 
relevant equation 

[12] 

is (5.112) of the book by Roberts* 
F^x,Q 2 ) = 2e 2 c asiQ \ +4ml) f dyg(y,Q 2 +4m 2 c )f(x/y,Q 2 ) (4a) 

^ Jax 

where 



f(z,Q 2 )=v Az 2 (l-z)-\, 



2m 



2ml 



_£ Z 2 (1 _ Z )j + ^ -z 2 (l-z) + -^z 2 (l - 3z) 



4m 



Q 4 



c z 3 



a = 1 + 4m 2 /Q 2 
v 2 = l-4m 2 c /[Q 2 (y/x-l)] 



(46) 



* Note that the formula Roberts gives for v should be for v 2 
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Figure 2: Charm cross section: pQCD calculation (thick lines) and the form (2) (thin lines). The 
data for Q 2 > are from ZEUS^L The photoproduction data are from Hit 9 ! and ZEUSt 10 !, who give 



references to the fixed-target data. The line in the lower right-hand corner of the photoproduction 
plot is our calculation for 6-quark production (for mass 4 GeV). 

For this leading-order calculation we again! 11 ] set Aqcd = 140 MeV. We have to choose the argument 
of a s and the scale \x of the gluon structure function; physical intuition leads us to take Q 2 + 4m 2 
for both, though it must be recognised that this is a mere guess. We need also to fix a value for m c . 
We find that 1.3 GeV gives good results: the thick lines in figure 2 are the output of the calculation, 
while the thin lines are the phenomenological fit (2). 



As we have said, our gluon distribution 
production we are modelling the streng" 
can be seen from the plots in figure 2, 
rise steeper than W 2e " at small W. To 
of the gluon distribution for values of x 
used (3) multiplied by (1 — x) n with n = 
cross section by less than 15% at W = 



is hard-pomeron dominated; when we use it to calculate charm 
th of the coupling of the hard pomeron to the charm quark. As 
the calculation also includes threshold effects which make the 
a small extent, these threshold effects depend on the behaviour 
that are beyond the small-x region where (3) is valid. We have 
: 5; changing n by one unit changes the charm photoproduction 
10 GeV and by 1% or less when W > 50 GeV. 
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Figure 3: Various gluon structure functions^] at Q 2 = 8 GeV 2 



At low Q 2 , and particularly for photoproduction, the magnitude of the cross section is sensitive to the 
value chosen for m c . Changing m c by 100 MeV away from our preferred value of 1.3 GeV changes 
the photoproduction cross section by more than 20% at the higher energies. We have not included 
any possible contribution from the hadronic structure of the photon, because its magnitude is so 
uncertain. Our calculations are consistent with it being small, but this may not be true^. If indeed 
it is small, one needs a structure function such as ours or the old to reproduce the steep 

Independence of the data^ at small Q 2 . The more modern MRST200lM and CTEQ6Mt 2 ] gluon 
structure functions are not large enough, and not steep enough, to reproduce the low-Q 2 data, as is 
obvious from figures 1 and 3. 

So, with /i 2 = Q 2 + 4m 2 and m c = 1.3 GeV, leading-order photn-gluon fusion gives a good description 
of the data for F 2 with hard-pomeron exchange alone, even down to Q = 0. At small enough x, the 
corresponding c-quark density is found to be almost identical with the hard-pomeron components of 
the densities of the light quarks. We previously® extracted the latter from the data for F 2 (x,Q 2 ) 
and! 11 ] showed that, for Q 2 greater than about 5 GeV 2 , they agree very well with 4-flavour zero-mass 
DGLAP evolution. It is 

standard[ 16 ][17][18] 

that photon-gluon fusion at small Q 2 must be matched 
to DGLAP evolution at large Q 2 . Our calculation has this property: the two agree over a large range 
of Q 2 values. 

We have not attempted to make a best fit to the charm-production data. In this, our approach is 
different from the so-called global fits 

[1,2]. 

Rather, our emphasis is on simplicity; we used only three 
free parameters in the fit to the small- x behaviour of F 2 (x,Q 2 ), and introduced no additional ones 
for F 2 (x,Q 2 ). The physics underlying our approach is very different from the conventional one and 
our successful application of pQCD to F 2 (x, Q 2 ) is a striking confirmation of the correctness of our 
extraction of the hard component from F 2 (x,Q 2 ). Most of F 2 (x,Q 2 ) at small x is the contribution 
from the soft component; this therefore plays the dominant role in the conventional approach, but it 
has not entered at all into the analysis we have described here. 

The first plot in figure 2 shows also our calculation for the 6-quark photoproduction cross section using 
rrib = 4 GeV. It is not inconsistent with a measurement of Hll^]. 

We now consider the proton's longitudinal structure function Fl(x, Q 2 ). We calculate this in leading- 
order pQCD. The relevant equation is (5.110) of the book by Roberts^]. However, we include the 
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effect of the mass m c in the contribution from c, c; this correction may be found in equation (E.3) of 
the review by Budnev et alt 2 ^ . So 

F L (x,Q 2 ) = G(x,Q 2 ) + ^^± r^(-) 2 F 2 ( y ,Q 2 ) (5) 

3vr J x y \yJ 

where the contribution of the charm quark to G(x, Q 2 ) is 

with a,v and L defined in (4b). We have again had to choose the argument of a s and the Q 2 -scale 
of the gluon structure function. We have made the same choice as before: Q 2 + Am 2 for each. Again 
this is a guess and the output is sensitive to it at low Q 2 . The light quarks contribute similarly, with 
m c replaced with 0. 

The HERA experiments measure the reduced cross section 

a r (x, y, Q 2 ) = F 2 (x, Q 2 ) - y _ y)2 F L (x, Q 2 ) (7) 

We have calculated this from our fit to F 2 and our gluon structure function (3); the results are shown 
in figure 4. The ranges of x and Q 2 shown are chosen because our fit to F 2 used only small- x data 
and because we foundt 11 ] that perturbative evolution only described it well for Q 2 greater than about 
5 GeV 2 , so our gluon distribution is not reliable for smaller values. In figure 5 we compare our 
calculated Fl with the values extracted by Hit 21 ! from their data. Separation of Fl from F 2 requires 
extrapolation and depends on some assumed parametrisation. 

In conclusion, in this paper we have continued our programme of reconciling the Regge and pQCD- 
evolution approaches to structure function data. We use two sets of data, the charm structure function 
and the longitudinal structure function. In our calculation of F% photon-gluon fusion at small Q 2 is 
matched to DGLAP evolution at large Q 2 . Our success in describing the charm and longitudinal 
structure functions provides further confirmation for the correctness of the two approaches and our 
understanding of how they fit together. We should warn, however, that in each case the extraction of 
the data from the raw measurements requires very large extrapolations. As MRST have observed M , it 
would be particularly useful to have good data for F L (x, Q 2 ), since this offers rather direct information 
about the gluon distribution. 

Our gluon structure function is larger and steeper at small values of x than is conventionally believed M t 2 L 
particularly at small Q 2 , but we have shown in this paper that there is experimental support for it. 
A less-steep gluon distribution will not explain the charm data at small Q 2 . Because our evolution 
procedure does not introduce spurious singulariies at N = into the splitting matrix, we conclude that 
the gluon distribution is larger than is usual in order to achieve the observed evolution of F 2 (x,Q 2 ). 

This research is supported in part by the EU Programme "Training and Mobility of Researchers", Network "Quantum 
Chromodynamics and the Deep Structure of Elementary Particles" (contract FMRX-CT98-0194), and by PPARC 
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Figure 4: Hit 21 ] and ZEUSt 22 ! data for the reduced cross section (7). The open squares are the values 
of F 2 (x,Q 2 ) extracted from these data by the two experiments. The data in (a) range from Q 2 = 5 
GeV 2 at the bottom to 60 at the top and in (b) from 6.5 to 120; in each case the data are separated 
by adding an extra 0.5 at successive values of Q 2 . At each Q 2 the upper line is our fitt 8 ] to F 2 (x,Q 2 ) 
(the ZEUS data in (b) were not available when the fit was made) and the lower line is the calculated 
reduced cross section. 
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Figure 5: HI data! 21 ! for the longitudinal structure function at various Q 2 values, with our pQCD 
calculations. The lower line on the Q 2 = 20 GeV 2 plot is the MRST2001 prediction M . 
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